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Abstract. Mirror matter is a stable self-collisional dark matter candidate. If exact mirror parity is a conserved symmetry of 
nature, there could exist a parallel hidden (mirror) sector of the Universe which has the same kind of particles and the same 
physical laws of our (visible) sector. The two sectors interact each other predominantly via gravity, therefore mirror matter 
is naturally "dark". Here I briefly review the cosmological signatures of mirror dark matter, as Big Bang nucleosynthesis, 
primordial structure formation and evolution, cosmic microwave background and large scale structure power spectra, together 
with its compatibility with the interpretation of the DAMA annual modulation signal in terms of photon-mirror-photon 
kinetic mixing. Summarizing the present status of research and comparing theoretical results with observations/experiments, 
it emerges that mirror matter is not just a viable, but a promising dark matter candidate. 
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INTRODUCTION 

Mirror matter is a stable self-interacting dark matter candidate that emerges if one, instead of (or in addition to) 
assuming a symmetry between bosons and fermions (supersymmetry), assumes that nature is parity symmetric. It is a 
matter of fact that the weak nuclear force is not parity symmetric. The main theoretical motivation for the mirror matter 
hypothesis is that it constitutes the simplest way to restore parity symmetry in the physical laws of nature. When Lee 
and Yang proposed the non-parity of weak interactions in 1956, they mentioned also the possibility to restore parity 
by doubling the number of particles in the Standard Model [1]. Thereby the Universe is divided into two sectors with 
the same particles and interactions, but opposite handedness, that interact mainly by gravity. In the minimal parity- 
symmetric extension of the Standard Model [2, 3], the group structure is G®G ', where G = SU(3) ®SU(2) <g)t/(l) and 
the prime (') denotes, as usual, the mirror sector. In this model the two sectors are described by the same lagrangians, 
but where ordinary particles have left-handed interactions, mirror particles have right-handed interactions. Besides 
gravity, mirror matter could interact with ordinary matter via the so-called kinetic mixing of gauge bosons, or via 
unknown fields that carry both ordinary and mirror charges. If such interactions exist they must be weak and are, 
therefore, negligible for many cosmological processes, but could be determinant for the detection of mirror dark 
matter in some contexts. Since photons do not interact with mirror baryons, or interact only via the weak kinetic 
mixing, mirror matter constitutes a natural candidate for the dark matter in the Universe. 

Many consequences of mirror matter for particle physics and astrophysics have been studied during the last decades. 
The reader can refer to J3l for a review of the history of mirror matter and a list of most relevant papers published 
before 2006. 

Like their ordinary counterparts, mirror baryons can form atoms, molecules and astrophysical objects such as 
planets, stars and globular clusters. However, even though the microphysics is the same in both sectors, the initial 
conditions, and then the macrophysics, of the mirror sector should be different, because the cosmology must be 
different. In particular, Big Bang nucleosynthesis (BBN) requires that the mirror sector has a lower temperature than 
the ordinary one This has implications for the thermodynamics of the early Universe and for the key 

cosmological epochs. Analytical results and numerical calculations of BBN SIM] show that the mirror sector should 
be helium dominated, and the abundance of heavy elements is expected to be higher than in the ordinary sector. 
Invisible stars made of mirror baryons are candidates for Massive Astrophysical Compact Halo Objects (MACHOs), 
which have been observed via microlensing events Mirror stars contain more helium and less hydrogen, 

and therefore have different properties than ordinary stars [12]. The accretion of mirror matter onto celestial objects 
Il3l[l^[l5tl and the presence of mirror matter inside compact stars, as for example neutron stars, could have interesting 
observable effects B16lll7ll . The consequences of mirror matter on primordial structure formation, cosmic microwave 
background (CMB) and large scale structure (LSS) distribution of matter have been studied HH [l9l l20l l2ll I22I l23ll . 
These studies provide stringent bounds on the mirror sector and prove that it is a viable candidate for dark matter. In 



addition, mirror matter provides one of the few potential explanations for the recent DAMA/LIBRA annual modulation 
signal I24ll25i[ 26 |l27ll . and possibly other experiments, as suggested for low energy electron recoil data from the CDMS 
collaboration 112811 . 

If the mirror (M) sector exists, then the Universe along with the ordinary (O) particles should contain their mirror 
partners, but their densities are not the same in both sectors. In fact, the BBN bound on the effective number of extra 
light neutrinos implies that the M sector has a temperature lower than the O one, that can be naturally achieved in 
certain inflationary models JH]. Then, two sectors have different initial conditions, they do not come into thermal 
equilibrium at later epoch and they evolve independently, separately conserving their entropies, and maintaining 
approximately constant the ratio among their temperatures. 

All the differences with respect to the ordinary world can be described in terms of only two free parameters: 

where T (T'), £2/, (Oj,), and s (s') are respectively the ordinary (mirror) photon temperature, cosmological baryon 
density, and entropy density. The lightest bounds on the mirror parameters are x < 0.7 and j3 > 1, the first one coming 
from the BBN limit and the second one from the hypothesis that a relevant fraction of dark matter is made of mirror 
baryons. 

As far as the mirror world is cooler than the ordinary one, x < 1, in the mirror world all key epochs (as are 
baryogenesis, nucleosynthesis, recombination, etc.) proceed in somewhat different conditions than in ordinary world. 
Namely, in the mirror world the relevant processes go out of equilibrium earlier than in the ordinary world, which has 
many far going implications. 

In this paper I consider the cosmological effects of the gravitational interaction between O and M matter, and I add 
the effects of photon-mirror-photon kinetic mixing only when showing the compatibility of the mirror dark matter 
interpretation of DAMA/LIBRA annual modulation signal with the present cosmological scenario. 



THERMODYNAMICS OF THE EARLY UNIVERSE 

In general, during most history of the Universe, we can approximate 
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where now q(T) and q'(T) are respectively the O and M entropic degrees of freedom, and the M photon temperature 
T'(T) is a function of the O one. This approximation is valid only if the temperatures of the two sectors are not too 
different, or otherwise if we are far enough from crucial epochs, like e + -e~ annihilation. When we investigate the 
range of temperatures interested by this phenomenon, as for example studying the BBN case, we need to exactly 
compute the trends of the thermodynamical quantities |7[]. 

Since ordinary and mirror sectors have the same microphysics, we may consider that the neutrino decoupling 
temperature Tq v is the same in both of them, that is Tqv = T^ v . We use this fact, together with the entropy conservation, 
to find equations which will give the mirror photon temperature T ! and the ordinary and mirror thermodynamical 
quantities corresponding to any values of the ordinary photon temperature T. From them it is possible to work out 
the total effective number of degrees of freedom (DOF) in both sectors, which can be, as common in the literature, 
expressed in terms of total effective number of neutrinos N v . The presence of the other sector, indeed, leads in both 
sectors to the same effects of having more particles. 

At temperatures we are interested in (below ~ 10 MeV) we can in general use the following equations derived from 
the conservations of separated entropies: 



22 _ \q e {T')+q y (T 1 



21 lq v \n 

22 lg e (T)+q r /T_ 
21 lq v \T V 



i 



3 



(3) 



TABLE 1. Effective N v in the ordinary sector. 



T(MeV) standard x = 0.1 x = 0.3 x = 0.5 x = 0.7 



5 3.00000 3.00063 3.04989 3.38430 4.47563 

0.005 3.00000 3.00074 3.05997 3.46270 4.77751 



TABLE 2. Effective N v in the mirror sector. 



T(MeV) 


x = 0.1 


x = 0.3 


x = 0.5 


x = 0.7 


5 


61432 


761.4 


101.3 


28.59 


0.005 


74011 


917.0 


121.4 


33.83 



and neglecting the entropy exchanges between the sectors (imposing x constant): 
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where s is the entropy density, a the scale factor, g, the entropic DOF of species i. 

At T ~ 7bv (Jp v ) ordinary (mirror) neutrinos decouple and soon after ordinary (mirror) electrons and positrons 
annihilate. 

The mirror world must be colder than the ordinary one and therefore the neutrino decoupling takes place before in 
the mirror sector. We can split the early Universe evolution in three phases. 

1) T > T Dv r. photons and neutrinos are in thermal equilibrium in both worlds, that is T v = T , Ty = T'. 

2) 7bv < T < T Dv i : at T ~ T Dv i M neutrinos decouple and soon after M electrons and positrons annihilate, raising the 
M photon temperature (T' ^ Ty). Nevertheless, O photons and neutrinos still have the same temperature (T = T v ). 

3) T < Tov- at T ~ Tov O neutrinos decouple and soon after O electrons and positrons annihilate, raising the O photon 
temperature (T ^T v ). 

Equations ([3]) and © are solved numerically in order to work out the total, ordinary and mirror numbers of entropic 
(q) and energetic (g) DOF at any temperature T. The corresponding number of neutrinos jV v is found assuming that 
all particles contributing to the Universe energy density, to the exclusion of electrons, positrons and photons, are 
neutrinos; in formula this means 
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In Table Q] I report the asymptotic numerical results before (T = 5 MeV) and after (T = 0.005 MeV) BBN for 
different values of x. I stress that the standard value A^ v = 3 is the same at any temperatures, while a distinctive feature 
of the mirror scenario is that the number of neutrinos raises with temperature and with x. Anyway, this effect is not 
a problem; on the contrary it may be useful since recent data fits give indications for a number of neutrinos at recent 
times higher than at BBN. 

It is possible to approximate the difference between the effective numbers of neutrinos before and after BBN process 
with the following expression: 



N V (T « T anne ±) -N V (T » T Dv ) =x 4 -^— 
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The effective number of neutrinos in the mirror sector can be worked out in a similar way, and the values (much 
higher than the ordinary ones) are reported in Table|2] 



BIG BANG NUCLEOSYNTHESIS 

As we have seen, the presence of the mirror sector can be parametrized in terms of extra DOF number or extra neutrino 
families; therefore, since the physical processes involved in BBN are not affected by the mirror sector, it is possible to 
use and modify a standard numerical code to work out the light elements production. 



TABLE 3. Light elements produced in the ordinary sector. 





standard 


x = 0.1 


x = 0.3 


x = 0.5 


x = 0.7 


4 He 


0.2483 


0.2483 


0.2491 


0.2538 


0.2675 


D/H (1(T 5 ) 


2.554 


2.555 


2.575 


2.709 


3.144 


3 He/H{W- 5 ) 


1.038 


1.038 


1.041 


1.058 


1.113 


7 Li/H (1CT 10 ) 


4.549 


4.548 


4.523 


4.356 


3.871 



TABLE 4. Light elements produced in the mirror sector (/3 = 5). 





x = 0.1 


x = 0.3 


x = 0.5 


x = 0.7 


4 He 


0.8051 


0.6351 


0.5035 


0.4077 


D/H (10~ 5 ) 


1.003- 10- 7 


4.838 -10- 4 


6.587-10- 3 


3.279- 10" 2 


3 He/H (10~ 5 ) 


0.3282 


0.3740 


0.4172 


0.4691 


7 Li/H(W- 10 ) 


1.996 - 10 3 


3.720- 10 2 


1.535 -10 2 


0.7962 ■ 10 2 



The number of DOF enters the program in terms of neutrino species number; this quantity is a free parameter, but 
instead of using the same number during the whole BBN process, I use the variable N v (T,x) numerically computed 
following the procedure described in the previous Section. The only parameter of the mirror sector which affects 
ordinary BBN is x; the baryonic ratio /3 does not induce any changes on the production of ordinary nuclides, but it 
plays a crucial role for the mirror nuclides production. 

In Table |3]I report the final abundances (mass fractions) of the light elements 4 He, D, 3 He and 7 Li produced in the 
ordinary sector at the end of BBN process (at T ~ 8 • 10~ 4 MeV) for several x values and compared with a standard 
scenario. We can easily infer that forx < 0.3 the light element abundances do not change more than a few percent, and 
the difference between the standard and x — 0.1 is of order 10~ 4 or less. 

Even mirror baryons undergo nucleosynthesis via the same physical processes than the ordinary ones, thus we 
can use the same numerical code also for the M nucleosynthesis. Mirror BBN is affected also by the second mirror 
parameter, that is the M baryon density (introduced in terms of the ratio j3 = Qi/Oj ~ 1 -f- 5), which raises the baryon 
to photon ratio rj' = j3x~ 3 rj. 

The results are reported in Table |4] which is the analogous of Table [3] but for a mirror sector with j3 = 5. We can 
see that BBN in the mirror sector is much more different from the standard than the BBN in the ordinary sector. 
This is a consequence of the high ordinary contribution to the number of total mirror DOF, which scales as ~ x~ 4 
(while in the ordinary sector the mirror contribution is almost insignificant, since it scales as ~ x 4 ). Hence, in this 
case the M helium abundance should be much larger than that of the O helium, and for x < 0.5 the M helium gives 
a dominant mass fraction of the dark matter of the Universe. This is a very interesting feature, because it means that 
mirror sector can be a helium dominated world, with important consequences on star formation and evolution 1 12], 
and other related astrophysical aspects. In particular, I recall the attention of the reader on the fact that the predicted 
dark matter composition, dominated by mirror helium, is exactly what is required by the proposed interpretation of the 
DAMA/LIBRA annual modulation signal in terms of mirror dark matter ll24ll25ll . 

Considering the different possible non-standard BBN scenarios, it would be of relevance a careful inspection to the 
effects of the presence of mirror particles on the predicted abundances of light elements. A more extended investigation 
in this sense is required in order to evaluate if it could help solving the still present problems on standard BBN due to 
the lithium anomalies. 



STRUCTURE FORMATION 

The important moments for the structure formation are related to the matter-radiation equality (MRE) and to the 
matter-radiation decoupling (MRD) epochs. The MRE occurs at the redshift 
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Therefore, for x <C 1 it is not altered by the additional relativistic component of the M sector. The mirror MRD 
temperature r d ' ec can be calculated following the same lines as in the O one, obtaining r d ' ec ss Td ec , and hence 

l+z dec ~x- 1 (l+z dec )^1100x- 1 , (8) 

so that the MRD in the M sector occurs earlier than in the O one. 

Moreover, for values x < x eq ~ 0.046 (i2,„/i 2 ) , the mirror photons would decouple yet during the radiation 
dominated period. This critical value plays an important role in our further considerations, where we distinguish 
between two cases: x > x eq and x < x eq . For typical values of Q. m h 2 we obtain x eq ~ 0.3. 

The relevant scale for gravitational instabilities is the mirror Jeans mass, defined as the minimum scale at which, 
in the matter dominated epoch, sub-horizon sized perturbations start to grow. In the case x > x eq (where the mirror 
decoupling happens after the matter-radiation equality) its maximum value is reached just before the M decoupling, 
and is expressed in terms of the O one as 

/ x 4 \3/2 

Kmax-r'^y^J -Mj, max , (9) 

which, for j3 > 1 and x < 1, means that the Jeans mass for the M baryons is lower than for the O ones, with implications 
for the structure formation. If, e.g., x = 0.6 and j3 = 2, then Mj ~ 0.03 Mj. We can also express the same quantity in 
terms of Q./,, x and j3, in the case that all the dark matter is in the form of M baryons, as 

MKfl dec )«3.2.1O 14 M ^- 1 / 2 (l+j3)- 3 / 2 (^ T ^- ? ) ' (£l h h 2 )- 2 ■ (10) 

For the case x < x eq , the mirror decoupling happens before the matter-radiation equality. In this case we obtain for the 
highest value of the Jeans mass just before decoupling the expression 

^(a d j«3.2ao i % ^-v2 (1+ ^ ) -3/ 2 ^y /2 ^_^y /2 (£2b/l 2 r 2 

In case x — x eq , the expressions (1 10b and (II Q , respectively valid for x > x eq and x < x eq , are coincident, as we expect. 
If we consider the differences between the highest mirror Jeans mass for the particular values x = x eq /2, x = x eq and 
x = 2x eq , we obtain the following relations: 

<max(-W2) « 0.005 Mj max (x eq ) ; 

<max(2*eq) « 64M J '. max (x eq ) . (12) 

Density perturbations in M baryons on scales M > Mj max , which enter the horizon at z ~ z e q, undergo uninterrupted 
linear growth. Perturbations on scales M < Mj max start instead to oscillate after they enter the horizon, thus delaying 
their growth till the epoch of M photon decoupling. 

As occurs for perturbations in the O baryonic sector, also the M baryon density fluctuations should undergo the 
strong collisional Silk damping around the time of M recombination, so that the smallest perturbations that survive the 
dissipation will have the mass 

M^[/(x)/2] 3 (j3i2,/2 2 )- 5 / 4 lO 12 M , (13) 

where f(x) — x 5 / 4 for x > x eq , and f(x) — (x/x eq ) 3//2 Xe q 4 for x < x eq . For x ~ x eq we interestingly obtain M' s w 
1Q 1 (D. b h 2 )- 5 M Q ~ 1O 1O M , a typical galaxy mass. 



COSMIC MICROWAVE BACKGROUND AND LARGE SCALE STRUCTURE 

In order to obtain quantitative predictions I computed numerically the evolution of scalar adiabatic perturbations in a 
flat Universe in which a significant fraction of mirror dark matter is present at the expenses of diminishing the cold 
dark matter (CDM) contribution and maintaining constant £l m . 

We have chosen a "reference cosmological model" with the following set of parameters: (Ot — Qi 7 h 2 — 0.023, £l m = 
0.25, Q.a = 0.75, n s = 0.97, h = 0.73. The dependence of the CMB and LSS power spectra on the parameters x and 




FIGURE 1. CMB (left) and LSS (right) power spectra for different values of x and co^, as compared with a reference standard 
model (solid line) and with observations. Models where dark matter is entirely due to M baryons (no CDM) are plotted in top panels 
for x = 0.3,0.5,0.7, while models with mixed CDM + M baryons (J3 = 1,2,3,4 ; x = 0.7) in bottom panels. 

j3 is shown in Fig. Q] The predicted CMB spectrum is quite strongly dependent on the value of x, and it becomes 
practically indistinguishable from the CDM case for x < x eq « 0.3. However, the effects on the CMB spectrum rather 
weakly depend on the fraction of mirror baryons. As a result of the oscillations in M baryons perturbation evolution, 
one observes oscillations in the LSS power spectrum; their position clearly depends on x, while their depth depends 
on the mirror baryonic density. Superimposed to oscillations one can see the cut-off in the power spectrum due to the 
aforementioned mirror Silk damping. 

In the same Figure our predictions can be compared with the observational data in order to obtain some general 
bound on the mirror parameters space. 

• The present LSS data are compatible with a scenario where all the dark matter is made of mirror baryons only if 
we consider enough small values of x: x < 0.3 « x eq . 

• High values of x, x > 0.6, can be excluded even for a relatively small amount of mirror baryons. In fact, we observe 
relevant effects on LSS and CMB power spectra down to values of M baryon density of the order Q/ b ~ Q.^,. 

• Intermediate values of x, 0.3 < x < 0.6, can be allowed if the M baryons are a subdominant component of dark 
matter, Q.b < £l' b < Qcdm- 

• For small values of x, x < 0.3, the M baryons and the CDM scenarios are indistinguishable as far as the CMB 
and the linear LSS power spectra are concerned. In this case, in fact, the mirror Jeans and Silk lengths, which 
mark region of the spectrum where the effects of mirror baryons are visible, decrease to very low values, which 
undergo non linear growth from relatively large redshift. 

Thus, with the current experimental accuracy, we can exclude only models with high x and high Q.' h . 



PHOTON-MIRROR-PHOTON KINETIC MIXING AND DAMA EXPERIMENT 



We have so far considered only gravitational interactions between O and M particles. Nevertheless, I mentioned in the 
Introduction the existence of other kind of weak interactions linking the two sectors. The most important of these is 
realized if O(M) particles couple to the M(O) particles via renormalizable photon-mirror-photon kinetic mixing 12911 . 
The mixing term in the Lagrangian is 

where F^ v = d^A v - d v A^ and F'^ v = d^A' v -d v A'^. This mixing enables 0(M) charged particles to couple to M(O) 
photons with charge ±ee. In the last years it has acquired more and more importance, since it has been shown that 
mirror dark matter is able to explain the positive dark matter signal obtained in the DAMA/Libra experiment [24], while 
also being consistent with the null results of the other direct detection experiments rf25ll . This interpretation in terms of 
mirror dark matter requires e ~ 10~ 9 , which is consistent with current laboratory and astrophysical constraints. The 
implications of such mixing for the early Universe have been recently studied in ref. B26I1 . where it has been shown 
in particular its consistency with constraints from ordinary BBN as well as more stringent constraints from CMB and 
LSS considerations. 

Assuming an effective initial condition T' <C T, we can estimate the evolution of T'/T in the early Universe, and also 
estimate the He' abundance. Photon-mirror-photon kinetic mixing can populate the mirror sector in the early Universe, 
via the process e + e~ — > e + 'e~'. The <? + ',<?~' will interact with each other via mirror weak and mirror electromagnetic 
interactions, populating the /, v' e) vL , v' z , and thermalizing to a common mirror sector temperature T'. 

Employing the usual expressions for the densities of relativistic fermions in thermal equilibrium and integrating 
with the initial condition T' = at T = Tj, we derive the evolution of the T'/T ratio as a decreasing function of T 



T -c 0.52 .10V/^-^V /4 . 0.164 e '/*(_L_V /4 . (15) 
T \ T Tj ) - g \T(MeV) J 



The last equivalence is obtained considering Tj lMeV and defining e_9 = e/10~ 9 . Clearly in all these computations 
we have neglected the change of the T'/T ratio due to the e + -e~ annihilation processes in both sectors (for a detailed 
study of this effect see ref. [7]), and the transfer from the mirror sector to the ordinary one (since T' < T). 

Note that the T'/T ratio freezes out when T < 2m e since the number density of e becomes Boltzmann suppressed 
and the process e + e~ — > e + 'e~' can no longer effectively heat the mirror sector. In addition, after e + '-e~' annihilation 
the effective number of mirror degrees of freedom decreases, and thus we estimate: 

T' 

— ~ 0.2 e l _'g for T < lMeV . (16) 

Constraints from ordinary BBN suggest that 8N V < 0.5 => T'/T < 0.6. A more stringent constraint arises from CMB 
and LSS, which suggest JH, H(J 111 T'/T < 0.3 and implies, from Eq. (16), that £ < 3 x 10~ 9 . This upper limit is 
compatible with the value required in order to explain the DAMA annual modulation signal 12511 . 

Given the above expression for the T'/T evolution, we can estimate the mirror helium mass fraction as a function of 
£. It has already been discussed in ref.(6l|8t] that, compared with the ordinary matter sector, we expect a larger mirror 
helium mass fraction if T' < T. Essentially, this is because the expansion rate of the Universe is faster at earlier times, 
which implies that the freeze out temperature of mirror weak interactions will be higher than that in the ordinary sector. 

The relevant time scales for BBN are the freeze out temperature of weak interactions Tw, and the "deuterium 
bottleneck" temperature T^ and time tpj. 

In the mirror sector, we have the same relations describing the BBN process, except that we change Tw — » T^ and 
fjV - ¥ i'n- Imposing the equality of weak interaction and Hubble rates Fw{T w ) = H(T/y), we obtain 

Tw^B-^e-^T^ 1 , (17) 
with B = 0.74 • 10 15 MeV and where we used Tw — ( -pr f, — ) . as computed from the equation Tw(Tw) — H(Tw). 

\ G' F Mpi J 

The quantity t' N can be estimated assuming T/j ~ TJv (since T^ depends only logarithmically on the baryon to photon 



ratio, that is larger in the M sector): 



f\ -4' 



-1/2 



M Pl _ R l/2 cT -l/2 t 

-~Y ~B el N t N , 

1 N 



(18) 



where we used t N ~ 0.3 g~ l/2 ^ and B- l e~ 2 T N > 1. 



Finally, using Eqs. ( TPTI i and ( 1181 , we can estimate the mirror helium mass fraction: 



2exp(-tL/x) 



2exp 



Y> ~7y'(t' t' ) ~ 1 v Nl ' 



_ B l/2 r -l/2 



Am B 2 / 7 



(19) 



Using the typical values T w ~ 0.8 MeV, 7]v ~ 0.07 MeV, Am ~ 1 .29 MeV, B ~ 0.74 • 10 15 MeV, x ~ 886.7 s, f w ~ 200 
s, we obtain: 



, 2exp[-2.3-10 7 e] 2exp[-2.3-10" 2 e_ 9 ] 
Ya ~ 1 + ^ [0.3 ■ 105 W] ~ l+exp r . 22 e% r 



(20) 



Thus, for e ~ 10 , we find F 4 ' ~ 0.87, that means a mirror sector largely dominated by helium. 



SUMMARY 

The current situation of the astrophysical research in presence of mirror dark matter is shown in Figure [21 with 
emphasis on the connections between the different theoretical studies and the predicted observable signatures. Solid 
lines mark what is already done (and is compatible with current observational and experimental constraints), while 
dashed ones mark what is still to do. For the last point a special importance is covered by star formation and N-body 
simulation studies, towards which people interested in mirror dark matter should address their efforts in the near future. 

At present, we have already investigated the early Universe (thermodynamics and Big Bang nucleosynthesis), and 
the process of structure formation in linear regime, that permit to obtain predictions, respectively, on the primordial 
elements abundances, and on the observed cosmic microwave background and large scale structure power spectra. In 
addition, we have studied the evolution of mirror dark stars, which, together with the mirror star formation, are nec- 
essary ingredients for the study and the numerical simulations of non linear structure formation, and of the formation 
and evolution of galaxies. Furthermore, in future studies they can provide predictions on the observed abundances of 
MACHOs and on the gravitational waves background. Using the results of mirror BBN and stellar evolution, we have 
verified that the interpretation, in terms of photon-mirror-photon kinetic mixing, of the DAMA annual modulation sig- 
nal is compatible with the inferred consequences on the physics of the early Universe. Ultimately, we will be able to 
obtain theoretical estimates of gravitational lensing, galactic dark matter distribution, and strange astrophysical events 
still unexplained (as for example dark galaxies, bullet galaxy, ...), that can be compared with observations. 

Concluding, the astrophysical and experimental tests so far used show that mirror matter is not just a viable, but a 
promising candidate for dark matter. However, we still need to obtain a complete picture of the Mirror Universe. 
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